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I.  INTRODUCTION 


A  major  goal  of  this  study  is  to  ascertain  the  ability  of  Molecular  Dynamics 
simulations  to  replicate  the  mechanical  properties  of  Parmax®  self-reinforced  polymers 
(SRPs).  These  polymers  consist  of  chains  of  para-  and  meta-  linked  phenylene  units 
with  benzoyl  group  substituents.  In  recent  years,  substituted  polyphenylenes  have 
emerged  as  a  new  class  of  compounds  with  excellent  mechanical  properties  (high 
modulus,  yield  strength  and  hardness).  They  also  display  excellent  thermal  stabilities 
and  processabilities  (low  Tg  and  soluble  in  many  organic  solvents). 

Prior  to  initiating  investigating  the  Parmax  SRPs,  we  studied  atactic 
Polypropylene  (aPP),  which  is  a  glassy  polymer  with  a  well  characterized  glass 
transition  temperature  and  reasonably  well  characterized  mechanical  properties. 


H2  H2  H2  Ha  H2  H2  H2  H2  H2 


We  used  our  simulations  on  this  polymer  to  develop  procedures  necessary  to 
determine  the  various  mechanical  moduli  as  well  as  the  yield  strength  in  amorphous 
polymeric  systems. 

Following  the  studies  on  aPP,  we  have  chosen  to  investigate  three  commercially 


available  benzoyl  substituted  polyphenylenes  (Ultem-1000,  Parmax-1200, 
Parmax-1500) 1  as  well  as  unsubstituted  poly(paraphenylene)  [PPP]  for  comparison. 


Ultem-1000  Parmax-1200  Parmax-1500  PPP 
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One  notes  from  the  structures  that, like  PPP,  Ultem-1000  has  100%  para  linkages  (as 
well  as  2,3-benzoyl  groups  on  all  monomers).  Parmax-1200  contains  only  85%  para 
linked  phenylenes  (and  85%  benzoyl  substitution),  and  Paramax-1500  is  only  50%  para 
linked  (and  50%  substitution). 

Below,  we  give  more  details  on  the  systems  investigated  (Sect.  II),  the 
computational  methods  used  in  the  simulations  (Sect.  Ill),  and  then  the  computed  vs. 
experimental  properties,  first  of  aPP  and  then  of  the  Parmax  SRPs  (Sect.  IV). 

II.  SYSTEMS  INVESTIGATED 

With  a  goal  towards  developing  the  protocols  required  to  compute  accurate 
elastic  moduli  and  yield  stresses  in  amorphous  polymers,  initial  investigations  in 
Summer  2004,  were  performed  on  the  well-studied  glassy  system,  atactic  polypropylene 
[aPP],  which  has  the  structure: 


Hj  h2  h2  h2  h2  h2  h2  h2  h2 


The  molecular  dynamics  cells  were  comprised  of  three  aPP  oligomers,  each  containing 
76  monomer  units  [689  total  atoms,  including  H3C-  and  -H  endcaps]. 

During  Summer  2005,  we  have  begun  analogous  investigations  of  three  Parmax 
Self-Reinforced  Polymers  (SRPs):  (A)  Ultem-1000  [U-1 000],  (B)  Parmax-1200 
[PX-1200],  (C)  Parmax-1500  [PX-1500],  in  addition  to  the  nonsubstituted  analogue  (D) 
Poly(paraphenylene)  [PPP], 
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The  MD  cells  were  designed  to  contain  a  minimum  of  >1500  atoms  and  have  cell 


lengths  of  at  least  25  A. 

U-1000:  3  strands 

PX-1200:  3  strands 

PX-1500:  4  strands 

PPP:  6  strands 


These  four  systems  were  comprised  of: 

25  monomers/strand  1 656  total  atoms 


25  monomers/strand 
25  monomers/strand 
25  monomers/strand 


1512  total  atoms 
1 608  total  atoms 
1512  total  atoms 


III.  COMPUTATIONAL  PROCEDURES 

Simulations  were  performed  using  the  Accelrys 2  Discover  molecular  dynamics 
program  and  associated  software  modules.  The  Accelrys  COMPASS® 3-6  force  field, 
designed  to  yield  accurate  values  for  both  gas  phase  and  liquid  phase  properties,  was 
used  for  all  the  simulations.  COMPASS  is  a  completely  atomistic  Class  II  force  field, 
and  the  potential  energy  function  contains  valence  terms  involving  bond  lengths,  angles, 
and  dihedrals  +  cross  terms.  As  with  other  modern  force-fields,  there  are  also 
out-of-plane  terms  designed  to  maintain  planarity  in  aromatic  systems,  including  the 
phenyl  rings  present  in  all  of  the  model  systems  investigated  in  this  work.  Other  than 
terms  in  the  potential  energy  function,  no  further  constraints  were  applied  to  the 
molecules  (i.e.  bond-lengths,  ring  geometries,  etc.,  were  not  frozen).  As  discussed  in 
the  original  references,3'6  the  nonbond  parameters  were  optimized  to  reproduce  molar 
volumes  and  enthalpies  of  vaporization  in  a  suite  of  test  molecules  under  ambient 
conditions  (usually  298  K  and  1  atm.  pressure),  in  which  aromatic  systems  were 
represented  by  benzene  and  toluene.5 
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Cubic  cells  containing  each  model  system  were  built  via  a  Monte  Carlo  algorithm 
as  implemented  in  the  Accelrys 1  Amorphous  Cell  module.  Following  minimization  to 
remove  close  contacts,  the  cells  were  heated  to  a  temperature  well  above  the  glass 
transition  temperature,  Tg,  and  were  equilbrated  for  XX  ns  using  NPT  dynamics. 
(Andersen  thermostat8  and  Berendsen  barostat9).  All  simulations  were  performed  with 
periodic  boundary  conditions  and  10  A  group-based  cutoffs.  The  step  size  was  1  fs. 

Equilibrated  cells  at  lower  temperatures  were  obtained  by  successive  cooling  at 
50  K  intervals  down  to  50  K,  and  equilibrating  using  NPT  dynamics  at  each 
temperature.  The  equilibrated  cell  volumes  from  these  runs  were  used  to  determine  the 
density  and  specific  volume  at  each  temperature. 

In  order  to  acquire  the  tabular  data  necessary  for  the  calculation  of  elastic  moduli 
(stresses  and  strains  in  each  cell  direction),  further  NPT  simulations  were  performed  at 
each  temperature  using  the  the  Parrinello-Rahman  method  (xx  ref),  in  which  an  NPT 
simulation  was  performed  with  Parrinello  pressure  control  and  XX  temperature  control. 
The  difference  between  this  run  and  earlier  NPT  simulations  is  that  the  x,y  and  z 
dimensions  of  the  box  are  permitted  independently,  allowing  the  calculation  of  strains 
(deformations)  along  the  three  cell  axes.  The  initial  box  dimensions  and  atomic 
coordinates  were  taken  from  the  final  frame  of  the  earlier  equilibration  runs.  The 
simulation  times  were  1100  ps,  and  the  first  100  ps  were  ignored  in  the  calculation  of 
elastic  properties  ( vide  infra). 

For  the  determination  of  stress  vs.  strain  curves,  a  BTCL  script  was  written  to 
perform  a  simulation  in  which  (a)  the  stress  in  the  x-direction  was  increased  (usually  in 
50  bar  increments)  up  to  a  maximum  of -2500  bar,  with  the  time  at  each  stress 
determined  by  the  stress-rate  (values  of  5,  3  and  1  bar/ps  were  used).  The  procedure 
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was  repeated  with  stresses  in  the  y  and  then  z  directions.  A  Matlab  script  was  then 
written  to  further  analyze  the  results  by  averaging  the  stress-strain  curves  in  the  three 
experiments. 

IV.  RESULTS  AND  DISCUSSION 
IVA.  Atactic  Polypropylene 

In  order  to  develop  and  test  the  procedures  for  determining  properties  of  the 
Parmax  SFPs,  we  chose  to  run  MD  simulations  on  a  simpler  polymer,  atactic 
polypropylene  (aPP).  This  system  (a)  is  known  to  be  completely  completely  amorphous 
(0%  crystallinity),  (b)  has  reasonably  well  characterized  volumetric  properties  (density, 
thermal  expansion  coefficients  and  glass  transition  temperature),10'12  and  (c)  its 
mechanical  properties  have  been  investigated  both  experimentally13  and  via  earlier  MD 
simulations.14  Two  independent  sets  of  simulations,  using  different  randomly  built 
amorphous  cells,  were  performed  on  aPP.  The  average  properties  (and  mean 
deviations)  are  presented  below. 

A.1  Volumetric  Properties  of  aPP 

Contained  in  Table  1A  are  the  available  experimental  specific  volumes  of  aPP11' 
13  as  well  as  volumes  (computed  from  NPT  runs)  for  one  of  the  simulations  over  a  range 
of  temperature  from  50  K  -  500  K.  Calculated  specific  volumes  from  the  two  simulations 
were  generally  very  close,  and  agreed  to  within  <0.5  %,  with  the  largest  differences 
<1.3%.  Calculated  and  experimental  volumes  are  also  displayed  in  Fig.  1  (filled  and 
open  symbols,  respectively).  One  observes  from  both  table  and  figure  that  there  is  very 
close  agreement  between  computed  and  experimental  values,  to  within  <0.7%  in  the 
high  temperature,  rubbery  regime.  The  authors  are  not  aware  of  any  experimental 
densities  in  the  glass  phase,  below  Tg. 
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Lines  tangent  to  the  computed  volumes  in  the  glass  and  rubbery  phases  were 
determined  from  the  four  lowest  and  highest  temperatures,  respectively,  and  are 
displayed  as  dotted  lines  in  Fig.  1  (displaced  downward  for  clarity  of  presentation). 

Using  the  standard  procedure,  the  glass  transition  temperature  was  determined  by  the 
intercept  of  the  two  lines.  From  the  data  in  the  figure,  it  was  found  for  the  first  of  the  two 
simulations  that  Tg(cal)  =  297  K;  the  second  run  (not  shown)  yielded  Tg(cal)  =  289  K. 
The  average  computed  glass  transition  temperature,  293  K,  is  16-18  K  above  two 
recent  determinations11,12  but  43  K  above  an  earlier  reported  value.  It  is  not  surprising 
that  the  computed  value  of  Tg  is  higher  than  experimental  values.  This  is  generally 
found  in  MD  simulations  as  a  result  of  the  extremely  high  effective  cooling  rates,  which 
elevate  the  computed  phase  transition  temperature. 

Shown  also  in  Table  1 A  and  Fig.  1  (dashed  curve)  are  specific  volumes 
calculated  from  from  QSPR  correlations,  Eqs.  [A.1]  and  [A.2],  with  required  parameters 
V(298)  taken  from  experiment  (Table  1A)  and  <xr(298)  determined  from  Eq.  [A.4]  (using 
Tg  =  275  K).  As  shown  in  the  figure  the  volumes  estimated  from  QSPR  agree  extremely 
well  with  both  experiment  and  the  MD  simulations  (to  within  at  high  temperatues,  and 
are  in  reasonable  qualitative  agreement  with  Vsp(MD)  in  the  low  temperature,  glassy 
regime  (deviations  average  1  %-2%  over  the  range  of  temperatures  in  the  glass). 
Interestingly,  the  glass  transition  temperature  predicted  from  the  intercept  of  lines  drawn 
tangent  to  the  QSPR  volumes  at  the  extremes  in  temperature  yield  a  value,  Tg  =  277  K, 
in  excellent  agreement  with  the  more  recent  experiments.11,12 

The  slopes  of  the  lines  tangent  to  the  computed  specific  volumes  (Fig.  1), 
(5Vsp/5T)p,  can  also  be  used  to  determine  the  thermal  expansion  coefficients  in  the 
glass  and  rubber  phases: 
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[1] 


1 


dV 


sp 


dT 


\ 

Jp 


Eq.  [1]  was  used  (with  average  volumes  over  the  temperature  range  spanned  by  each 
tangent  line)  to  determine  expansion  coefficients,  aG  and  ocr,  in  the  two  phases.  The 
average  values  (and  mean  deviations)  for  the  two  simulations  are  given  in  Table  1 B, 
together  with  available  experimental  results.  It  may  be  that  aG(MD)  lies  very  close  to 
the  middle  of  the  range  of  expansion  coefficients  published  for  aPP.  The  agreement 
between  o.r(MD)  and  aR(Exp)  is  also  good.  The  computed  expansion  coefficient  lies 
within  5%-10%  of  the  three  experimental  values  in  the  high  temperature  phase. 

Thermal  expansion  coefficients  computed  from  tangent  lines  to  the  QSPR 
volumes  (not  shown  in  the  figure)  are  contained  in  the  last  column  of  Table  1  A. 
Alternatively,  one  can  use  the  empirical  correlations  between  a(T)  and  Tg  (Eqs.  [A.3] 
and  [A.4]).  For  temperatures  representing  the  average  of  the  low  and  high  temperature 
regimes,  values  of  aG  and  aR  (125  K  and  425  K,  respectively),  obtained  from  the 
equations  are  in  excellent  agreement  with  those  computed  from  the  equations  and  from 
slopes  of  the  tangent  lines.  The  difference  between  QSPR  and  MD  values  of  aR  are 
insignificant  (differing  by  <  5%).  There  is  a  more  substantial  difference  between 
aG(MD)  and  aG(QSPR)  [.  Unfortunately,  because  of  the  large  range  in  experimental 
determinations,11'13  the  relative  accuracy  of  the  two  estimation  methods  cannot  be 
ascertained. 
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A.2  Elastic  Moduli  of  aPP 


The  six  elements  of  the  strain  tensor  (<7j)  are  related  to  the  six  stress 
components  (sj)  by  the  "stiffness  matrix"  (Cy)  via  the  transformation:15,16 
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Because  the  stiffness  matrix  is  symmetric  (i.e.  Cjj  =  Cy),  there  can  be  as  many  as  21 
independent  elements.  However,  the  symmetry  present  in  most  systems  will  reduce  the 
number  of  independent  terms  dramatically.  Of  specific  interest  to  us  is  the  matrix  for 
amorphous,  isotropic  systems,  in  which  case  all  matrix  elements  can  be  cast  in  terms  of 
the  two  Lame  coefficients,  X  and  p,  and  is  of  the  form:15,16 
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In  turn,  the  tensile  modulus  (E),  bulk  modulus  (B),  shear  modulus  (G)  and  Poisson's 
ratio  (v)  are  functions  of  the  Lame  coefficients,  and  given  by: 
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Probably  the  most  versatile  method  for  utilizing  a  Molecular  Dynamics  simulation  to 
determine  elements  of  the  stiffness  matrix  and,  hence,  the  elastic  moduli,  is  the 
"fluctuation  method,"  introduced  by  Parrinello  and  Rahman.17  They  proved  the 
individual  elements  of  C  are  related  to  correlations  in  strain  tensor  elements  via  the 


relation: 


where  T  and  T  and  <V>  are  the  temperature  (K)  and  averaged  cell  volume.  <£j£j> 


[8] 


represents  the  average  of  strain  tensor  components  oyer  all  frames  in  the  simulation's 
production  run.  More  recently,  Gusev  et  a/.18  pointed  out  that  the  original  fluctuation 
expression  for  Cy  is  slowly  convergent,  oftentimes  requiring  very  lengthy  simulations. 
They  proposed  the  the  following  alternative  formula  for  Cy,  which  converges  more 


rapidly : 


[9] 


To  obtain  the  requisite  stress  and  strains,  we  performed  ~1  ns  NPT  (1,000,000 
frames)  simulations  using  the  Parrinello-Rahman  barostat 17  and  Andersen  thermostat.8 
The  stiffness  matrix  elements  were  obtained  from  the  appropriate  averages,  using  Eqs. 
[8]  and  [9].  A  typical  matrix,  at  100  K  using  the  strain-strain  correlation  formula  is: 
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A  comparison  of  the  numerical  data  with  the  matrix  elements  expected  for  an 
isotropic  sample  (Eq.  [3])  reveals  that  the  results  are  in  good  qualitative  agreement;  e.g. 
Cii  «  C22  ~  C33,  C12 «  C13 «  C23,  and  the  remaining  elements  are  significantly  lower. 

One  obtains  E,  B,  G  and  v  via  a  least-squares  fit  of  the  non-zero  elements  to  obtain  X 
and  p,  followed  by  application  of  Eqs.  [4]-[7]. 

The  temperature  dependence  of  the  elastic  constants  computed  from  the  MD 
simulations  are  displayed  in  Table  2;  E,  B  and  G  are  also  plotted  in  Fig.  2.  Values 
presented  represent  those  obtained  from  the  strain-strain  correlation,  Eq.  [8].  Although 
not  shown,  constants  from  the  stress-strain  correlation  (Eq.  [9])  are  very  close  in 
magnitude.  One  observes  from  both  table  and  figure  that,  all  three  moduli  decrease 
monotonically  with  rising  temperature.  This  is  expected  and  reflects  softening  of  the 
material  at  higher  temperatures.  In  addition,  the  MD  results  for  the  elastic  moduli  vary 
in  the  order:  B  >  E  >  G,  as  is  found  experimentally  for  most  elastic  materials.  In 
addition,  as  seen  from  the  last  column  of  Table  2,  values  of  the  Poisson  ratio  lie 
between  0.35  and  0.40,  and  vary  only  modestly  with  temperature.  This  is  common 
behavior  for  amorphous  polymers  well  below  Tg. 

The  author  is  not  aware  of  any  mechanical  experiments  that  have  been 
performed  on  atactic  polypropylene.  In  one  of  the  common  references  commonly 
quoted  for  polypropylene  mechanical  properties  in  compilations,19'21  the  results  are 
actually  for  the  isotactic  stereoisomer.22  A  fairly  comprehensive  study  of  the 
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temperature  dependence  of  tensile  moduli  of  amorphous  polypropylene  has  been 
reported  by  Sauer,  et  a/.13  For  one  of  their  samples  (PP-4  in  the  reference),  they  report 
0%  crystallinity  (from  IR/X-ray  results),  but  do  not  indicate  the  tacticity  of  the  sample. 
However,  to  obtain  an  estimate  of  the  qualitative  accuracy  of  our  results,  we  have 
extracted  values  of  E(T)  [via  numerical  digitization  of  data  in  Fig.  2  of  their  paper],  which 
are  displayed  in  Table  2  [E(exp)]  and  plotted  in  Fig.  2  [open  circles].  It  was  found  that 
the  qualitative  agreement  between  E(MD)  and  E(exp)  is  fairly  good;  values  from  MD  are 
typically  15-30%  lower  than  experiment.  The  temperature  dependence  of  the  computed 
tensile  moduli  [dE(MD)/dT  =  -9.1  MPa  K'1]  is  also  close  the  the  experimental  value 
[dE(exp)/dT  =  -10.4  MPa  K'1].  The  level  of  agreement  between  calculated  and 
measured  moduli  is  particularly  satisfying  given  uncertainties  in  the  tacticity  of  the 
sample,  the  much  greater  molecular  weight  of  the  experimental  specimen  [Mw(exp)  « 
3x1 05  g  mol'1  vs.  MW(MD)  »  3,200  g  mol'1],  and  the  well-known  dependence  of 
measured  mechanical  properties  upon  the  method  of  sample  preparation  and 
processing  history.  Theodorou  and  Suter 14  have  also  reported  the  results  of  an 
atomistic  molecular  mechanics  simulation  of  the  elastic  moduli  of  aPP.  One  cannot 
assign  a  temperature  to  the  results  obtained  in  that  study  because,  although  the  model 
was  static  (i.e.  at  0  K),  the  amorphous  cell  was  built  at  a  density  corresponding  to  233 
K.  Nevertheless,  a  qualitative  of  the  two  sets  of  results  is  interesting.  As  shown  in  the 
last  row  of  Table  2,  their  computed  values  for  E,  B  and  G  all  lie  in  the  same  range  as 
our  results  between  the  temperatures  of  100  K  and  150  K,  from  which  we  may  conclude 
that  the  results  of  the  two  investigations  are  in  very  reasonable  qualitative  agreement. 

Finally,  it  is  instructive  to  compare  elastic  moduli  predicted  from  QSPR  with 
values  computed  from  our  MD  simulations.  These  results,  determined  from  Eqs. 
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[A.6]-[A.8],  are  plotted  in  Fig.  2  also  plotted  in  Fig.  2.  One  may  see  from  the  figure  that 
both  the  magnitudes  and  temperature  dependence  of  the  QSPR  agree  quite  well  with 
values  computed  from  our  MD  simulations.  Specifically,  B(QSPR)  agrees  with  B(MD)  to 
within  -1  %  to  -21  %,  with  similar  levels  of  agreement  for  E  (-6%  to  +1 1  %)  and  G  (-4%  to 
+  13%). 

A.3  Yield  Stress  of  aPP 

As  discussed  in  an  earlier  section,  computer  experiments  were  performed  in 
which  the  the  uniaxial  stress  along  a  given  (x,  y,  or  z)  axis  in  the  cell  was  increased  at  a 
constant  rate  and  the  resulting  strain  was  measured.  This  method  of  studying 
stress-strain  relationships  in  polymers  was  first  developed  by  Brown  and  Clarke,23  and 
their  or  similar  methods  have  been  employed  successfully  in  a  number  of  MD 
investigations  of  the  yield  behavior  in  solid  polymers  [e.g.  Refs. 24'29].  Displayed  in  Fig. 
3A  are  the  resultant  stress  vs.  strain  curves  for  aPP  at  three  temperatures  when  the 
stress  was  increased  at  a  rate  of  5  bar  ps'1;  similar  curves  were  found  for  lower  rates  of 
stress  increase.  One  observes  the  general  phenomenon  in  plastic  solids  that  at  low 
stresses,  the  strain  increases  linearly  with  applied  “true  stress”30  (elastic  behavior). 
However,  at  higher  stresses,  the  sample  begins  to  exhibit  a  very  large  increase  in  strain 
(i.e.  catastrohic  failure)  with  little  or  no  further  increase  in  the  stress  levels. 

Because  the  specific  shapes  of  stress  vs.  strain  curves  vary  considerably 
between  samples  (and  even  for  the  same  sample  at  different  samples),  there  are  a 
variety  of  definitions  of  the  yield  stress  (oy)  of  a  solid,15  and  of  the  procedures  by  which 
Gy  is  derived  from  a  specific  curve.  Qualitatively,  though,  one  may  view  the  yield  stress 
as  that  value  beyond  which  the  plastic  solid  begins  to  exhibit  liquid-like  flow  behavior 
(i.e.the  point  at  which  the  stress-strain  curve  approaches  a  horizontal  increase  in 
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strain).31  In  order  to  obtain  semi-quantitative  estimates  of  the  yield  stress  (ay)  and  yield 
strain  (ey)  in  our  aPP  samples,  we  utilized  a  common  procedure15  in  which  tangents 
were  drawn  in  both  the  elastic  (low  s)  and  yield  (high  s)  region  of  the  curves;  sy  and  ay 
were  determined  from  the  intersection  of  the  two  lines.  Because  the  high  yield  line  is 
close  to  horizontal,  values  of  ay  are  expected  to  be  more  accurate  than  sy  (this  is  true  in 
all  graphical  procedures  used  to  calculate  the  yield  stress). 

Numerical  estimates  of  the  yield  stress  and  strain  at  3  temperatures  and  3  rates 
of  stress  increase  are  contained  in  Table  3.  The  values  represent  the  average  of  two 
independent  runs  (using  independently  built  and  equilibrated  cells),  with  the  mean 
deviation  between  runs  in  parentheses. 

As  seen  clearly  in  both  Fig.  3A  and  in  Table  3,  for  a  given  rate  of  stress  increase, 
the  yield  stress  decreases  markedly  with  increasing  temperature.  This  is  expected 
intuitively,  because  of  the  “softening”  of  the  plastic  at  higher  temperature,  and  has  been 
observed  in  other  MD  investigations  of  yield  behavior.23'29 

One  also  expects  the  stress-strain  curves  and  resultant  yield  stresses  to  be  a 
function  of  the  rate  of  stress  (or  strain)  increase  on  the  sample.  We  have  plotted  the 
curves  for  three  rates  of  stress  increase  at  a  given  temperature  (100  K)  in  Fig.  3B. 
Indeed,  one  finds  that  the  yield  stress  is  lowered  as  the  rate  of  increase  is  decreased; 
this  is  also  found  in  Table  3  (although  the  values  of  ay  at  200  K  are  approximately  equal 
for  5  bar  ps'1  and  3  bar  ps'1).  Again,  this  observation  is  expected  on  a  physical  basis 
because  decreased  rates  of  stress  increase  provide  the  time  required  for  the  polymer 
segments  to  reequilibrate  in  response  to  increased  stress,  thus  promoting  additional 
yielding  of  the  sample.  The  same  trend  has  been  observed  in  other  MD  simulations  of 
polymer  stress-strain  behavior.23'29 
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To  the  author’s  knowledge,  the  only  reported  yield  stress’s  of  polypropylene  have 
been  on  the  isotactic  polymer  or  on  samples  of  unknown  tacticity.10  Furthermore,  it  is 
well  known  that  measured  values  ay  on  the  same  polymer  can  vary  markedly  and  are 
dependent  upon  quantities  such  as  molecular  weight,  sample  preparation  method  and 
processing  history.  However,  it  is  still  possible  to  obtain  an  estimate  of  the  quality  of 
the  simulation  results  by  referring  to  the  well-known  empirical  correlation  that  the  yield 
stress  in  a  large  number  of  polymeric  materials  is  proportional  to  the  tensile  modulus: 

Oy(T)«K.E(t)  [11] 

with  reported  values  of  the  empirical  proportionality  constant  constant,  k,  ranging  from 
0.02520  to  0.028. 21 

Displayed  in  Fig.  4  are  plots  of  the  yield  stress  in  aPP  as  a  function  of 
temperature  for  each  of  the  three  stress  rates  (5,  3  and  1  bar  ps_1).  For  comparison,  a 
plot  of  the  experimental  tensile  moduli  (Eexp)13  over  a  similar  temerature  range  is 
superposed  (with  the  axis  for  Eexp  on  the  right  side  of  the  figure).32  It  is  very  satisfying  to 
note  that  the  slopes  of  the  computed  yield  stresses  are  qualitatively  quite  similar  to  that 
of  the  experimental  tensile  moduli,  indicating  that  the  MD  simulations  have  correctly 
captured  the  expected  temperature  dependence  of  the  yield  stresses  in  aPP. 

Quantitatively,  we  find  that  the  ratios,  ay(MD)  /  Eexp,  average  about  0.050-0.060, 
which  is  approximately  twice  the  ratio  expected  from  the  empirical  correlation  found  for 
other  polymers  (Eq.  11);  i.e.  the  simulated  values  of  ay  are  too  high  by  approximately  a 
factor  of  two.  In  actuality,  this  result  is  not  at  all  surprising  because,  by  necessity,  the 
rate  of  stress  increase  in  the  simulations  are  many  orders  of  magnitude  greater  than 
those  used  experimentally.  Thus,  one  expects  the  computed  values  of  ay  to  be  higher 
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than  experiment;  this  has  also  been  observed  in  other  MD  simulations  of  yield  stresses 
in  glassy  polymers. 23'29 

Unlike  many  properties,  no  accurate  semi-quantitative  relationships  have  been 
developed  to  predict  the  dependence  of  yield  stresses  on  stress  (or  strain)  rate. 
However,  it  has  been  reported  that,  theoretically,  one  expects  that  ay  should  be  an 
increasing  linear  function  of  the  logarithm  of  the  strain-rate  increase  (sR):19'21 


—  a  -Y  f5  •  log  10  )  [12] 

with  (3  >  0.  As  discussed  in  the  computational  section,  the  method  used  in  determining 
the  stress-strain  behavior  utilizes  the  rate  of  increase  of  the  stress  (not  strain)  as  an 
independent  variable.  However,  by  assuming  an  approximately  linear  plot  of  stress  vs. 
strain  (Fig.  3)  up  to  the  yield  point  (ey,  ay),  we  can  estimate  the  strain  rates  (sr)  from  the 
stress  rates  (SR)  via  following  relation: 

[13] 


£  R  ~ 


sy  ■  SR  1012 


where  the  factor  1012  converts  SR  to  bar  s'1,  and  the  units  of  £r  are  %»s'1.  As  noted 
above,  values  of  ey  are  subject  to  considerable  error.  Therefore,  in  applying  Eq.  12,  we 
have  used  the  average  value  (over  the  three  stress  rates)  of  sy.  Values  of  sr  derived  in 
this  manner  are  given  in  the  last  column  of  Table  3. 

Shown  in  Fig.  5  are  plots  of  ay  vs.  log(sR)  at  all  three  temperatures  (100  K,  200  K, 
300  K).  One  can  observe  from  the  figure  that,  to  within  the  scatter  generated  by 
inherent  uncertainties  in  determination  of  ay,  the  yield  stresses  do  increase 
approximately  linearly  with  logarithm  of  the  strain  rate,  as  predicted  by  theory. 

In  summary,  we  have  found  that  computed  yield  stresses  decrease  with  rising 
temperature,  and  increase  at  higher  strain  rates.  Both  trends  are  expected  intuitively 
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and  observed  experimentally.  Furthermore,  the  computed  values  of  oy  are  within  a 
factor  of  two  of  values  estimated  from  an  empirical  correlation  with  the  tensile  modulus, 
and  exhibit  the  same  temperature  dependence  as  E,  as  predicted  by  the  correlation. 
IVB.  Parmax  Self-Reinforced  Polymers 

As  noted  in  the  introduction,  the  structures  of  the  Parmax  SRPs  and  of  PPP  are 
given  by: 


85%  50% 


Ultem-1000  Parmax-1200  Parmax-1500  PPP 

The  order  of  increasing  backbone  rigidity  (i.e.  percent  of  para  linkages)  is: 

Parmax-1500  <  Parmax-1200  <  Ultem-1000  =  PPP.  The  order  of  increasing  percent  of 
benzoyl  substitution  is:  PPP  <  Parmax-1500  <  Parmax-1200  <  Ultem-1000. 

Below,  we  first  discuss  volumetric  properties  [Sect.  IVB.1],  then  mechanical 
moduli  [Sect.  IVB.2]  and  then  yield  stresses  [Sect.  IVB.3].  Where  appropriate,  results 
on  the  Parmax  SRPs  are  compared  with  experiments  and  with  trends  in  molecular 
structure. 

B.1  Volumetric  Properties  of  Parmax  SRPs 

From  NPT  equilibration  simulations,  values  for  the  cell  dimensions  were 
determined,  which  yielded  densities  and,  hence,  specific  volumes  (vsp  =  1/p)  for  each 
polymer  as  a  function  of  temperature.  A  representative  plot  of  specific  volume  versus 
temperature  is  given  for  Parmax-1200  in  Fig.  6.  The  computed  specific  volume  at  room 
temperature  is  in  reasonable  agreement  with  experiment  (within  1.5%)  at  the  one 
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temperature  reported.  Values  of  the  glass  transition  temperature  were  computed  in  the 
standard  manner  as  the  intersection  of  the  low  and  high  temperature  tangent  lines  in  the 
specific  volume  plots.  Calculated  values  of  Tg  are  given  in  Table  4,  together  with 
experimental  values  for  the  three  Parmax  SRPs.33  It  may  be  seen  that  the  three 
commercial  polymers  exhibit  very  little  variation  of  Tg(exp)  with  structure  (only  a  total  of 
a  10  K  range).  One  finds  too  that  Tg(MD)  is  higher  than  Tg(exp)  by  24  to  57  K.  This  is 
not  surprising  and  has  been  found  elsewhere.  The  higher  values  of  Tg(MD)  is  attributed 
directly  to  the  much  higher  rates  of  sample  cooling  in  the  simulations  (by  many  orders  of 
magnitude)  than  in  the  experimental  measurements. 

There  is  no  obvious  trend  in  the  computed  glass  transition  temperatures  with 
molecular  structure.  For  example,  the  two  highest  values  of  Tg(MD)  are  for  PPP  and 
Parmax-100,  which  have  the  highest  and  lowest  number  of  para  linkages,  and  thus,  the 
lowest  and  highest  chain  flexibility,  respectively.  Too,  the  unsubstituted  PPP  sample 
exhibits  a  calculated  transition  temperature  within  6  K  of  the  value  for  Ultem-1000, 
which  has  100%  benzoyl  substitution. 

B.2  Elastic  Moduli  of  Parmax  SRPs 

Tabulated  in  Table  5  and  displayed  in  Fig.  7  are  (a)  the  experimental 
compressive  moduli  [in  MPa]  and  (b)  experimental  compressive  strengths  [in  bar]  in 
the  three  commercial  Parmax  SRPs  at  several  temperatures.  We  note  that  the 
compressive  modulus  should,  in  principle,  be  approximately  equal  the  tensile  modulus 
(E),  whereas  the  compressive  strength  is  generally  greater  than  the  tensile  yield 
strength  (by  perhaps  20-50%). 15’19'21 

From  Table  5A  and  Fig.  7A,  one  finds  that  the  moduli  vary  in  the  order: 
Ultem-1000  <  Parmax-1500  <  Parmax-1200.  One  sees  the  same  order  in  the 
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compressive  strengths  (Table  5A  and  Fig.  7B);  however,  from  the  figure,  the  strengths 
of  the  three  polymers  are  close  to  equal  at  a  temperature  just  below  400  K.  This 
temperature  is  approximately  30-40  K  below  Tg  in  the  systems.  However,  it  is  not  clear 
why  the  compressive  strengths  of  three  polymers  of  differing  structure  should  approach 
each  other  at  the  glass  transition  temperature. 

It  is  of  interest  that  the  order  of  increasing  experimental  compressive  moduli  and 
strengths  do  not  correspond  to  any  apparent  trend  in  structure.  For  example,  the  order 
of  increasing  molecular  backbone  rigidity,  as  indicated  by  percent  para  linkages  would 
be  Parmax-1500  (50%  para)  <  Parmax-1200  (85%  para)  <  Ultem-1000  (100%  para). 
Note  that  the  trend  in  increasing  benzoyl  substitution  is  the  same.  Thus,  the  trend  in 
these  structural  properties  does  not  correspond  to  the  experimental  order  of  mechanical 
properties,  in  which  the  Parmax-1500  has  properties  midway  between  those  of 
Parmax-1200  and  Ultem-1000. 

We  have  used  the  fluctuation  strain-strain  method  [Eq.  8  above]  to  determine  the 
mechanical  Tensile  (E),  Bulk  (B)  and  Shear  (G)  moduli  +  Poisson’s  ratio  (v)  in  PPP  and 
the  three  Parmax  SRPs.  The  results  are  displayed  in  Table  6  and  the  moduli  are 
plotted  in  Fig.  8.  As  shown  in  Table  6D,  Poisson’s  ratio  (v)  generally  lies  in  the  range 
0.35-0.45  expected  for  amorphous  polymers  at  temperatures  well  below  Tg. 19,21 
Generally,  the  dependence  of  v  on  temperature  is  very  slight  (until  one  reaches  close  to 
Tg)  and,  indeed,  the  values  of  this  parameter  appear  to  remain  approximately 
independent  of  T  to  within  scatter  in  the  computed  results. 

From  Table  6  and  Fig.  8A,  it  is  clear  that  the  tensile  moduli  (E)  of  all  four 
polymers  decreases  with  increasing  temperature.  This  is  the  same  trend  observed  for 
aPP  ( vide  supra)  and  expected  of  the  mechanical  moduli  in  amorphous  solids.15,19,21  A 
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similar  decrease  in  the  shear  modulus  (G)  [Fig.  8B]  and  bulk  modulus  (B)  [Fig.  8C]  is 
observed  for  all  four  systems  at  higher  temperatures. 

One  can  see  from  Tables  6A-6C  that,  for  a  given  polymer,  the  mechanical  moduli 
vary  in  the  order:  B(T)  >  E(T)  >  G(T).  As  noted  in  the  earlier  section  on  aPP,  this  is  the 
trend  generally  found  in  glassy  polymers  below  the  glass  transition  temperature. 

Inspection  of  Fig.  8A  refeals  that  the  calculated  tensile  moduli  vary  in  the  order: 
E(Ultem-IOOO)  >  E(Parmax-1200)  >  E(Parmax-1 500)  >  E(PPP)  at  most  temperatures. 
The  same  order  is  found  for  the  shear  moduli  (G)  of  the  four  polymers  [Fig.  8B].  On  the 
other  hand,  there  is  no  clear  cut  trend  in  computed  bulk  moduli  (B)  in  the  four  systems 
[Fig.  8C].  In  retrospect,  this  is  perhaps  not  surprising  because  B,  as  a  measure  of 
volume  decrease  with  increasing  isotropic  compression,  is  inherently  different  from  E 
and  G.both  of  which  are  measures  of  increased  elongational  strain  upon  the  application 
of  linear  stresses.  One  also  notes  that  the  observed  order  of  decreasing  tensile  and 
shear  moduli  is  the  same  as  the  trend  in  decreasing  benzoyl  substitution:  Ultem-1000 
[100%]  >  Parmax-1200  [85%]  >  Parmax-1500  [50%]  >  PPP  [0%].  While  this 
observation  is  interesting,  it  is  unfortunate  that  because  the  relative  backbone  rigidities 
of  the  polymers  (as  measured  by  fractional  meta  linkage)  varies  irregularly  in  the  series, 
one  may  not  conclude  that  the  trend  of  decreasing  mechanical  moduli  is  a  direct  results 
of  the  amount  of  side-chain  substitution. 

As  noted  earlier,  it  is  a  reasonable  approximation  to  assume  that,  for  small 
stresses,  the  compressive  and  tensile  moduli  are  close  to  equal. 15,19,21  Therefore,  we 
can  compare  our  results  for  the  tensile  moduli  with  experimental  compressive  moduli  for 
the  three  Parmax  SRPs.  This  is  shown  in  Fig.  9  (experimental  data  are  represented  by 
open  symbols). 
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One  first  notices  that,  in  general,  experimental  compressive  moduli  are 
significantly  greater  than  computed  tensile  moduli  in  Parmax-1200  and  Parmax-1500, 
by  factors  of  2  to  3,  although  experimental  and  calculated  values  for  Ultem-1000  are 
rather  close.  One  sees,  further  that  among  the  SRPs,  experimental  compressive  moduli 
vary  in  the  order:  Parmax-1200  >  Parmax-1500  >  Ultem-1000.  Thus,  although  the  the 
order  of  predicted  moduli  from  the  MD  simulations  is  correct  for  Parmax-1200  and 
Parmax-1500,  the  computed  moduli  of  Ultem-1000  are  predicted  to  be  the  highest 
among  the  SRPs,  which  is  precisely  the  opposite  from  the  trend  found  experimentally. 

As  discussed  earlier,  the  various  SRPs  differ  both  in  their  percent  substitution 
and  percent  of  backbone  flexibility.  Therefore,  it  is  not  possible  at  this  time  to  ascertain 
whether  the  results  indicate  that  there  are  errors  in  the  predictions  for  the  rigid  polymer 
(Ultem-1000),  which  lies  closest  to  experiment,  or  from  overestimation  of  the 
mechanical  moduli  for  systems  with  flexible  backbones  (Parmax-1200  and 
Parmax-1500).  In  the  last  section,  we  discuss  proposed  computational  experiments  to 
study  the  effects  of  benzoyl  substitution  and  backbone  rigidity  separately. 

B.3  Yield  Stress  of  Parmax  SRPs 

Uniaxial  stress  experiments  (see  Sect.  Ill:  Computational  Procedures)  were 
conducted  on  the  three  Parmax  SRPs  and  PPP  at  temperatures  ranging  from  50  K  to 
400  K  using  a  rate  of  stress  increase  of  5  bar  ps'1  in  order  to  determine  the  stress-strain 
behavior  in  these  systems.  In  the  earlier  section  on  yield  behavior  in  aPP,  it  was 
discussed  that  values  for  the  yield  stress,  oy  were  determined  from  intersection  of  the 
linear  regions  of  the  curves  for  low  and  high  stress  behavior  (see  Fig.  3).  However,  it 
was  found  that  the  stress-strain  curves  in  the  Parmax  SRPs  did  not  exhibit  the  same 
regular  behavior;  curves  for  the  same  polymer  at  different  behaviors  and  for  different 
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polymers  had  different  shapes.  Therefore,  we  employed  an  alternative  procedures  to 
estimate  the  yield  stress.  We  measured  cy  as  the  value  of  the  stress  at  which  the  strain 
reached  a  constant  percent.  Of  various  values  used,  the  most  regularly  trends  were 
found  when  values  of  ay  were  determined  as  the  stress  when  strain  reached  10%. 

Values  of  ay  estimated  in  this  manner  are  given  in  Table  7  for  the  four  systems 
as  a  function  of  temperature;  they  are  also  plotted  in  Fig.  10  (solid  symbols).  There 
were  some  points  inconsistent  with  the  surrounding  data.  These  are  given  in  bold  italics 
and  the  table  and  are  not  plotted. 

It  is  satisfying  to  see  from  the  figure  that,  in  general,  the  yield  stress  decreases 
as  the  temperature  rises.  This  was  found  in  aPP  and  is  expected  because  of  polymer 
“softening”  as  one  approaches  the  glass  transition  temperature. 

One  also  observes  that,  although  there  are  exceptions,  for  a  given  temperature 
the  yield  stresses  of  the  various  polymers  vary  in  the  order: 

cy(U  Item- 1000)  >  ay(Parmax-1200)  >  ay(Parmax-1 500)  >  ay(PPP).  This  precisely  the 
same  trend  observed  for  the  tensile  (and  shear)  moduli  (Fig.  8A)  and,  as  noted  above, 
is  the  same  as  the  order  of  decreasing  benzoyl  substitution. 

In  the  earlier  section  on  yield  stress  in  aPP,  it  was  discussed  that  there  is  an 
empirical  linear  relationship  between  yield  stresses  and  tensile  moduli  in  polymeric 
solids  [Eq.  (11)].  We  have  computed  values  for  the  ratio,  K=cry(T)/E(T),  in  the  SRPs 
(being  sure  to  convert  the  two  quantities  to  the  same  units).  The  values  of  k  range  from 
-0.05  to  0.09  (with  a  few  values  being  a  bit  higher),  which  is  approximately  2-3  times 
the  experimentally  observed  ratios.21  However,  it  is  satisfying  to  note  that  there  is  a 
clear  correlation  between  computed  values  of  these  two  mechanical  properties. 
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It  is  of  interest  to  compare  calculated  yield  stresses  in  the  Parmax  SRPs  with 
experimental  compressive  strengths  in  these  systems,  keeping  in  mind  that  (a) 
compressive  strengths  are  typically  20-50%  higher  than  yield  stresses  and  (b)  we  have 
only  rough  estimates  of  ay  in  these  systems  due  to  the  comparatively  irregular  behavior 
of  the  stress-strain  curves.  The  experimental  compressive  strengths  are  also  plotted  in 
Fig.  10  (open  symbols).  One  sees  that  the  calculated  and  measured  quantities  are  in 
satisfying  qualitative  agreement  with  each  other.  However,  it  can  be  seen  that,  whereas 
ay(U Item- 1000)  calculated  from  MD  is  highest  of  the  SRPs,  this  polymer  has  the  lowest 
compressive  strength  experimentally.  This  discrepancy  is  precisely  the  same  that  was 
found  when  comparing  computed  tensile  and  experimental  compressive  moduli  in  the 
last  section.  Once  again,  because  the  percent  of  benzoyl  substitution  and  backbone 
rigidity  are  both  varied  independently  in  these  systems,  it  is  not  possible  to  ascertain  the 
basis  of  the  deviation  between  theory  and  experiment,  at  this  time. 

V.  Further  Studies 

The  investigation  on  the  volumetric  properties,  elastic  moduli  and  yield  behavior 
in  atactic  polypropylene  is  close  to  complete,  and  the  data  yield  a  coherent  picture  on 
the  mechanical  properties  of  this  system.  On  the  other  hand,  there  are  many  questions 
which  remain  about  the  elastic  and  yield  properties  of  the  Parmax  SRPs.  A  number  of 
studies  are  planned  which  can  address  these  issues  and  yield  a  better  understanding  of 
the  capability  of  MD  simulations  to  model  the  mechanical  behavior  in  these  systems. 
Specifically,  we  plan  studies  of  the  following  type. 

1.  Simulations  using  independently  built  cells.  All  results  reported  in  this  study 
were  obtained  by  measurements  of  a  single  cell  of  each  polymer.  When  one  is 
computing  properties  in  a  polymeric  systems  containing  only  3-4  strands  of 
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polymer,  then  it  is  possible  that  the  results  will  differ  using  different  cells  with 
varying  structure.  The  experiments  should  be  repeated  with  at  least  one 
additional  cell  (and  preferably  two  or  more)  to  see  if  the  computed  mechanical 
properties  are  cell  independent. 

2.  Dependence  on  cell  size.  In  our  earlier  studies  of  fluid  properties,  we  found  that 
cells  with  lengths  of  ~25  A  were  adequate  to  replicate  macroscopic  behavior 
(using  PBCs,  of  course).  However,  this  may  not  be  true  in  polymeric  solids. 
Therefore,  we  will  perform  selected  experiments  with  larger  cells  (35-40  A  cell 
length)  to  determine  whether  the  results  remain  the  same. 

3.  Stress  Rate  dependence  of  yield  stresses.  It  was  noted  in  Sect.  IVB.3  that  the 
stress  strain  curves  of  the  Parmax  SRPs  yielded  irregularly  shaped  curves, 
unlike  earlier  results  on  aPP.  We  believe  that  this  may  be  a  result  of  the 
commparatively  high  rate  of  increase  of  the  applied  stress  (5  bar  ps'1). 

Therefore,  experiments  on  these  systems  will  be  repeated  at.  lower  rates  of 
stress  increase  of  3  and  1  bar  ps'1  to  determine  if  the  stress-strain  curves  exhibit 
more  regular  behavior  under  the  milder  conditions. 

4.  Comparison  with  predictions  from  QSPR.  It  was  found  in  our  pilot  studies  on 
aPP  that  empirical  Quantitative  Structure-Property  Relations  yielded  predicted 
volumetric  and  mechanical  properties  that  were  in  reasonably  good  agreement 
with  those  derived  from  the  MD  simulations.  Because  QSPR  furnishes  a  very 
computationally  inexpensive  method  of  property  prediction,  we  will  determine  if 
these  correlations  perform  adequately  in  predicting  the  various  properties  in  the 
Parmax  SRPs  as  well. 
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5.  Conformational  Properties.  It  is  well  known  that  polymer  solids  distort  under 
stress.  We  plan  to  determine  conformational  properties  including  end-to-end 
distances,  radii  of  gyrations,  dihedral  angle  distributions  and  transition  rates, 
orientational  distribution  parameters  and  atomic  pair  correlation  functions  in  the 
glassy,  rubbery  and  fluid  regimes  to  monitor  changes  in  these  parameters  in  the 
varioius  phases.  In  particular,  the  conformational  properties  in  the  glass  will  be 
measured  under  varying  stress  conditions  to  determine  the  correlation  between 
these  measures  of  molecular  conformation  and  the  computed  mechanical 
properties.  We  are  particularly  interested  in  comparing  dihedral  angle 
distributions  and  transition  rates  in  the  various  SRPs  to  determine  the  effect  of 
the  presumed  backbone  flexibility  (from  different  percent  meta  linkages)  on  these 
parameters. 

6.  Independent  variation  of  backbone  flexibility  and  benzoyl  substitution.  In 

the  above  discussion  of  the  elastic  moduli  and  stress  behavior  in  the  Parmax 
SRPs,  we  noted  that  an  unambiguous  determination  of  the  effect  of  structural 
features  on  the  mechanical  properties  was  not  possible  because  both  the 
fractional  substitution  and  chain  flexibilities  were  varied  concurrently  in  the 
commercial  polymers  (undoubtedly  as  a  result  of  synthetic  procedures).  With 
MD  simulations,  we  can  vary  these  structural  features  independently.  For 
example  we  will  determine  the  effects  of  benzoyl  substitution  on  mechanical 
properties  by  calculations  on  PPP  polymers  (100%  para  linkages)  with  varying 
percentages  of  benzoyl  group  substitution.  Similarly,  experiments  on 
unsubstituted  polyphenylenes  with  varying  percents  of  meta  linkages  will  enable 
us  to  determine  how  the  backbone  flexibility  affects  the  computed  properties. 


24 


Experiments  of  this  type  will  enable  us  to  assess  the  effects  of  substitution  and 
flexibility  independently. 

Appendix:  QSPR  estimation  of  polymer  properties 

The  application  of  Quantitative  Structure-Property  Relationships  (QSPRs)  to 
predict  physical  properties  in  polymeric  systems  has  been  the  subject  of  much  research, 
(Refs.  19,  20,  21  and  references  therein)  and  parameterized  semi-empirical  equations 
have  been  developed  which  permit  the  estimation  of  various  volumetric,  thermal, 
mechanical  and  other  properties  (e.g.  optical,  electrical,  solution,  etc.)  with  a  reasonable 
degree  of  accuracy.  We  have  utilized  several  of  these  derived  equations  to  predict  the 
values  of  various  properties  which  have  been  computed  from  the  MD  simulations. 
Volumetric  Properties 

Bicerano21  has  presented  a  series  of  equations  to  predict  volumes  (specific  or 
molar)  in  polymers.  For  systems  in  which  Tg  <  298  K,  one  has  the  relations: 


V(T)  =  0.15^(298) 


T-T 


1.427;  +44.7 


-+ 


V (298)[l  +  aR  (298)(r^  -  298)] 


T  <  T 


[A.  1  ] 


V(T)  =  K (298 )[l  +  «^(298)(r-298)]  T>  Tg  [A.2] 

V(298)  and  oir(298)  are  the  room  temperature  volume  and  thermal  expansion 
coefficient,  respectively. 

As  shown  by  Seitz,20  the  thermal  expansion  coefficients  in  the  glass  and  rubber 
can  be  predicted  by  the  relations: 
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Elastic  Properties 

Seitz20  has  developed  correlations  for  the  bulk  modulus  (B)  and  Poisson's  Ratio 
(v)  in  amorphous,  glassy  polymers  (we  use  the  form  of  the  equations  presented  in 
Bicerano's  book.21  The  temperature  dependence  of  v  is  given  by: 


r  A  rji 

v(T)  =  v0  +  —  {O.OOOI  63  +  exp[o.459(r  —  2^  —13)]} 
Tg 


[A.  5] 


vo  is  the  Poisson's  ratio  at  T  =  0  K,  and  can  be  determined  if  one  has  an  experimental 
(or  estimated)  value  of  v(T)  at  some  higher  temperature. 

The  temperature  dependence  of  the  bulk  modulus  can  be  predicted  from  the 
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expression. 


B(T)  =  S.2333Ecohl 


5^/(0)4 
.  Vmol(T)5 
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[A.  6] 


In  this  equation,  VmOi(0)  and  Vm0i(T)  are  the  molar  volumes  at  T  =  0  K  and  the 
temperature  of  interest,  and  ECOh  i  is  the  Fedors  estimate  (via  group  contributions)  of  the 
cohesive  energy.34  As  written,  if  Vm0i  and  ECOh  i  are  given  in  cm3  mol"1  and  J  mol'1, 
respectively,  then  the  units  on  B(T)  are  MPa. 

In  amorphous  glassy  solids,  there  are  only  two  independent  elastic  constants. 
Therefore,  if  v(T)  and  B(T)  are  known  experiment  or,  in  this  case,  estimated  from 
Eqs.  [A.5]  and  [A.6],  the  tensile  and  shear  moduli  (E  and  G)  can  be  obtained  from  the 
relations:15,35 
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E(T)  =  3B(T)[l-2v(T)] 

[A.7] 

G(T)  =  3  1  ~  2v(T^  b(T) 

2  1  +  v(T) 

[A.8] 

It  should  be  noted  that  the  dependence  of  Poisson's  ratio  on  temperature  is  very  modest 
until  one  approaches  temperatures  close  to  Tg  (e.g.  T  >  Tg-  20  K).  Therefore,  from  Eqs. 
[A.7]  and  [A.8],  one  expects  an  approximately  constant  proportionality  between  the 
three  moduli,  B,  E  and  G. 

Yield  Stress 

There  have  been  several  correlations  developed  between  the  ultimate  strength  of 
a  glassy  polymer,  as  measured  by  the  yield  stress  (ay)  and  the  elastic  moduli.  It  has 
been  found  that,  at  low  strain  rates,  there  is  a  direct  proportionality  between  ay(T)  and 
E(T):20,21 

<ry(T)  =  K-E(T)  [A.  9] 

Two  commonly  used  values  of  the  proportionality  constant  are  k  =  0.25  (Ref.  20)  or 
k  =  0.28  (Ref.  21). 

It  has  also  been  found  that  when  the  rate  of  strain  increase  on  a  polymeric  solid 
is  increased,  the  yield  stress  rises.  A  very  approximate  proportionality  is  given  by:21 


cr  (T)  =  k-E(T)  +  P- In 
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where  Vs  and  e  are  empirical  constants.  Thus,  one  expects  a  proportionality  between 


cry  and  ln(s)  at  high  strain  rates. 
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Table  1.  Volumetric  Properties  of  aPP. 


A.  Specific  Volumes 


T 

[K] 

Vsp[exp] 

[cm3  g*1] 

vj(An 

[cm3  g*1] 

Vsp[QSPR] 

[cm3  g'1] 

50 

1.092 

1.068 

100 

1.108 

1.088 

150 

1.125 

1.108 

200 

1.141 

1.129 

250 

1.163 

1.149 

298 

1.1 7b 

300 

1.175° 

1.184 

1.178 

350 

1.213° 

1.216 

1.218 

400 

1.256° 

1.264 

1.258 

450 

1.301 

1.298 

463 

500 

1.31 d 

1.346 

1.338 
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Table  1.  (Cont’d.) 


B.  Thermal  Expansion  Coefficients 
and  Glass  Transition  Temperatures 


Quantity 

Exp. 

MDxe 

QSPR 

104aG  [tC1] 

1.9a 

2.81 

3.67 

2.7f 

(0.13) 

[3.66]h 

3.9g 

104aR  [K*1] 

6.9b 

6.53 

6.26 

7.1d 

(0.13) 

[6.30]‘ 

Ui 

CM 

H 

CO 

3 

2509 

293 

277 

275b 

(4) 

277d 

Footnotes 

a)  Values  from  one  of  two  MD  runs. 

b)  Ref.  12. 

c)  Ref.  13.  Numerical  values  computed  from  numerical 
digitization  of  data  in  Fig.  3  of  reference. 

d)  Ref.  11. 

e)  Results  are  the  average  of  two  simulations. 

values  in  parentheses  are  the  mean  deviation  from  the  avers 

f)  Refs.  14. 

g)  Ref.  10. 

i)  Calculated  at  125  K  from  Eq.  [A.3],  Assumed  Tg  =  275  K. 
i)  Calculated  at  425  K  from  Eq.  [A.4].  Assumed  Tg  =  275  K. 
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Table  2.  Elastic  Moduli  of  aPP. 


T 

[K] 

E(exp)a 

[MPa] 

E(MD)b 

[MPa] 

B(MD)b 

[MPa] 

G(MD)b 

[MPa] 

v(MD)b 

100 

3870 

3350 

4320 

1220 

0.37 

(100) 

(10) 

(40) 

(<0.01 ) 

150 

3510 

2540 

3340 

920 

0.38 

(230) 

(10) 

(90) 

(0.01) 

200 

3040 

2310 

2880 

850 

0.37 

(300) 

(220) 

(120) 

(0.01) 

250 

2510 

1700 

2180 

620 

0.39 

(260) 

(10) 

(110) 

(0.02) 

C 

2890 

3320 

1065 

0.36 

Footnotes 

a)  Ref.  13.  Values  obtained  by  numerical  digitization  of  data  in  Fig.  3  of  reference. 

b)  Values  represent  average  of  the  two  simulations,  using  the  Strain-Strain 
correlation  equation  (Eq.  [8]). 

c)  Ref.  14.  Moduli  were  calculated  in  absence  of  kinetic  energy  (i.e.  at  0  K),  but  using 
cells  built  to  the  experimental  density  at  233  K. 

Moduli  represent  the  average  value  obtained  from  two  methods  in  Table  1 
of  the  reference. 
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Table  3.  Temperature  and  Stress  Rate  Dependence  of 
Yield  Stress  (cry)  and  Yield  Strain  (ey)  in  aPP.a 


T  Stress  Rate 

[K]  [bar  ps'1] 

C7y 

[bar] 

8y 

[%] 

£r 

[%  s-1] 

100  5 

2340 

7.6 

1.62E+10 

(3) 

(0.8) 

3 

2100 

7.7 

1.10E+10 

(70) 

(0.3) 

1 

2030 

7.3 

3.61  E+09 

(10) 

(0.6) 

200  5 

1770 

9.2 

2.61  E+10 

(40) 

(0.1) 

3 

1800 

10.6 

1.77E+10 

(30) 

(1.8) 

1 

1390 

8.9 

6.41  E+09 

(80) 

(2.5) 

300  5 

1400 

9.8 

3.49E+10 

(30) 

(0.6) 

3 

1300 

9.4 

2.16E+10 

(20) 

(3.2) 

1 

990 

7.4 

7.45E+09 

(10) 

(0.8) 

Values  represent  the  average  of  two  simulations. 

Quantities  in  parentheses  are  the  mean  deviation  from  the  average. 
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Table  4.  Experimental  and  Computed  Glass  Transition 
Temperatures  in  Parmax  SRPs. 

Species  Tg(exp)  Tg(MD) 

_ _  [K] _  [K] 

PPP  -  483 

U  Item-1 000  438  489 

Parmax-1200  433  457 

Parmax-1500  443  500 


a)  Ref.  33. 

b)  Mississippi  Polymer  Technologies,  Inc. 
Data  provided  to  AFRL/ML. 
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Table  5.  Experimental  Mechanical  Properties  of  Parmax  SRPs.a 
A.  Compressive  Modulus 


T 

Ultem-1000 

Parmax-1200 

Parmax-1500 

[K] 

[MPa] 

[MPa] 

[MPa] 

219 

4620 

9500 

7720 

298 

3790 

8840 

6810 

394 

3190 

6590 

5500 

B.  Compressive  Strength 

T 

Ultem-1000 

Parmax-1200 

Parmax-1500 

[K] 

[bar] 

[bar] 

[bar] 

219 

1860 

2730 

2510 

298 

1520 

2240 

1900 

394 

980 

1050 

1030 

a)  Missippii  Polymer  Technologies,  Inc. 
Data  provided  to  AFRL/ML 
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Table  6.  Computed  Elastic  Constants  of  Parmax  SRPs.a 
A.  Tensile  Modulus  (E) 


T 

PPP 

Ultem-1000 

Parmax-1200 

Parmax-1500 

[K] 

[MPa] 

[MPa] 

[MPa] 

[MPa] 

100 

2438 

4427 

3277 

3458 

200 

1742 

3799 

4348 

2657 

300 

1162 

2777 

2348 

1907 

400 

652 

2558 

2412 

1101 

B.  Bulk  Modulus  (B) 

T 

PPP 

Ultem-1000 

Parmax-1200 

Parmax-1500 

[K] 

[MPa] 

[MPa] 

[MPa] 

[MPa] 

100 

7693 

5366 

6255 

5127 

200 

5455 

4996 

5284 

4873 

300 

3887 

2956 

4234 

4273 

400 

1995 

3290 

3048 

2992 
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Table  6.  (Cont'd.) 


C.  Shear  Modulus  (G) 


T 

[K] 

PPP 

[MPa] 

Ultem-1000 

[MPa] 

Parmax-1200 

[MPa] 

Parm  ax-1 500 

[MPa] 

100 

842 

1624 

1160 

1246 

200 

602 

1383 

1595 

943 

300 

401 

1034 

834 

669 

400 

226 

933 

882 

383 

D.  Poisson's  Ratio  (v) 

T 

PPP 

Ultem-1000 

Parmax-1200 

Parmax-1500 

[K] 

100 

0.447 

0.363 

0.413 

0.388 

200 

0.447 

0.373 

0.363 

0.409 

300 

0.450 

0.343 

0.408 

0.426 

400 

0.446 

0.370 

0.368 

0.439 

a)  Computed  using  the  strain-strain  correlation  (Eq.  8). 
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Table  7.  Calculated  Yield  Stress  at  10%  Strain  in  Parmax  SRPs.a,b 

T 

[K] 

PPP 

[bar] 

Ultem-1000 

[bar] 

Parm  ax-1 200 

[bar] 

Parmax-1500 

[bar] 

50 

1770 

1898 

2868 

1796 

100 

1385 

2604 

2780 

1472 

150 

1219 

-- 

—  ' 

2196 

200 

1160 

2348 

— 

1563 

250 

— 

2318 

2127 

1440 

300 

1038 

1918 

1809 

1155 

350 

1095 

2129 

1264 

1457 

400 

760 

1910 

1763 

1104 

a)  Entries  with  a  indicate  that  the  run  "crashed"  before  the  strain  reached  10 

b)  Entries  in  bold  italics  are  inconsistent  with  surrounding  data.  These 
points  were  not  plotted  in  the  figure. 
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Figure  Captions 


1.  Specific  Volume  of  aPP. 

•  -  Calculated  from  MD  simulation  [A  represents  one  of  two  runs] 

O  -  Ref.  R2  (Eck) 

□  -  Ref.  R5  (Sauer) 

A  -  Ref.  R3  (Dlubek) 

- Calculated  from  QSPR. 

. -  Tangents  to  MD  volumes  in  glass  and  rubber  regions  [lines  have  been 

displaced  downward  by  0.04  cm3  g'1  for  clarity]. 

2.  Elastic  Moduli  of  aPP. 

■  -  B(MD);  Bulk  modulus  computed  from  MD  simulation. 

•  -  G(MD);  Shear  modulus  computed  from  MD  simulation. 

•  -  E(MD);  Tensile  modulus  computed  from  MD  simulation. 

O  -  E(exp);  Experimental  tensile  modulus;  Ref.  13 

- -  B(QSPR);  Bulk  modulus  from  QSPR. 

. -  G(QSPR);  Shear  modulus  from  QSPR. 

. E(QSPR);  Tensile  modulus  from  QSPR. 

3.  Stress  vs.  Strain  Curves  in  aPP. 

A.  Dependence  on  Temperature  (Stress  Rate  =  5  bar  ps'1). 

B.  Dependence  on  Stress  Rate  (T  =  100  K). 

4.  Temperature  Dependence  of  Yield  Stress  in  aPP. 

O - -  CTy  (SR  =  5  bar  ps'1) 

□  . <jy  (SR  =  3  bar  ps'1) 

O . -  CTy  (SR  =  1  bar  ps'1) 

•  -  E(exp) 

5.  Dependence  of  Yield  Stress  on  log(Strain  Rate)  in  aPP 

•  - -  T  =  100  K 

■  - T  =  200  K 

•  . -  T  =  300  K 

The  lines  through  the  points  were  drawn  qualitatively. 

6.  Specific  Volume  of  Parmax-1200. 

•  -  Calculated  from  MD  simulation 

O  -  Experimental  (Bars  represent  error  limits) 

. -  Tangents  to  MD  volumes  in  glass  and  rubber 
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7.  Experimental  Mechanical  Properties  of  Parmax  SRPs 

(A)  Experimental  Compressive  Moduli 

(B)  Experimental  Compressive  Strengths 

•  -  U  Item- 1000 

■  -  Parmax-1200 

•  -  Parmax- 1500 

8.  Calculated  Mechanical  Properties  of  Parmax  SRPs 

(A)  Calculated  Tensile  Moduli 

(B)  Calculated  Shear  Moduli 

(C)  Calculated  Bulk  Moduli 

A - -  PPP  [MD] 

•  - Ultem-1000  [MD] 

■  . -  Parmax-1200  [MD] 

Parmax-1500  [MD] 

9.  Comparison  of  Calculated  and  Experimental  Moduli  in  Parmax  SRPs 

•- . Ultem-1000  [MD] 

■  . -  Parmax-1200  [MD] 

Parmax-1500  [MD] 

O  -  Ultem-1000  [exp] 

□  -  Parmax-1200  [exp] 

O  -  Parmax-1500  [exp] 

10.  Calculated  Yield  Stresses  (at  10%  Strain)  of  Parmax  SRPs 

A - -  PPP  [MD] 

•  - Ultem-1000  [MD] 

■  -  Parmax-1200  [MD] 

♦----  Parmax-1500  [MD] 

O  -  Ultem-1000  [exp] 

□  -  Parmax-1200  [exp] 

O  -  Parmax-1500  [exp] 
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Fig.  3B:  Stress  vs.  Strain  Curve  in  aPP 
Dependence  on  Stress  Rate 
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Fig.  4:  Temperature  Dependence  of  Yield  Stress 


Fig.  5:  Dependence  of  Yield  Stress 
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Tensile  Modulus  (E)  [MPa] 


Fig.  8A:  Calc'd.  Tensi 
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Shear  Modulus  (G)  [MPa] 


Fig.  8B:  Calc'd.  Shear  Moduli  of  Parmax  SRPs 
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Fig.  8C:  Calc'd.  Bulk  Moduli  of  Parmax  SRPs 
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Fig.  10:  Calc'd.  Yield  Stresses  (at  10%  Strain) 
of  Parmax  SRPs 
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